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Ablation  Modeling  for  Dynamic  Simulation  of  Reentry 

Vehicles 

David  Doman*  and  William  Blake^ 

Air  Force  Research  Laboratory /VAC A,  Wright- Patters  on  AFB,  OH  45433-7531 


I.  Introduction 

The  collection  of  methods  described  in  this  manuscript  can  be  used  in  a  dynamic  simulation  to  provide 
estimates  of  the  mass  properties  and  aerodynamic  forces  and  moments  as  a  reentry  vehicle  ablates  due  to 
aerodynamic  heating  .  Dynamic  simulations  of  reentry  vehicles  are  required  to  perform  trajectory  optimiza¬ 
tion  and  to  analyze  the  robustness  of  guidance  and  control  systems.  Reentry  vehicle  trajectories  vary  widely 
depending  on  mission  requirements  and  vehicle  design.  Vehicles  that  experience  relatively  low  peak  heating 
with  reusable  thermal  protection  systems  such  as  the  shuttle,  experience  little  if  any  ablation.  At  the  other 
extreme,  ballistic  reentry  vehicles  and  interplanetary  probes  can  experience  very  high  peak  heat  loads  that 
cause  the  thermal  protection  material  to  ablate.  A  number  of  vehicle  characteristics  change  as  a  result  of  ab¬ 
lation.  The  mass  properties  of  the  vehicle  change  due  to  the  loss  of  material  and  the  aerodynamic  forces  and 
moments  acting  on  the  vehicle  change  as  a  result  of  the  an  ablating  outer  mold  line  (OML).  These  changes 
can  affect  aerodynamic  as  well  as  guidance  and  control  system  performance.  A  number  of  techniques  exist 
that  allow  one  to  measure  the  changes  to  the  OML1,2  during  experimental  testing.  Empirical  methods  are 
described  in  this  paper  that  can  be  used  to  translate  limited  test  data  into  a  rough,  but  representative  model 
that  can  be  used  to  estimate  the  effects  of  ablation  on  a  vehicle’s  ability  to  follow  a  prescribed  trajectory 
and  on  guidance  and  control  performance  and  robustness. 

Section  II  proposes  a  method  to  estimate  the  location  of  points  on  the  OML  that  lie  outside  of  the  vehicle 
plane  of  symmetry  where  test  measurements  are  typically  obtained.  These  points  can  be  used  to  generate  a  3 
dimensional  grid  that  is  suitable  for  use  with  aerodynamic  prediction  codes.  The  out-of-plane  grid  points  are 
estimated  using  a  radial  interpolation  method  that  can  accommodate  asymmetries  that  occur  between  the 
windward  and  leeward  surfaces  of  the  vehicle  .  Formulae  for  the  geometric  properties  of  radially  interpolated 
ablated  nose  cone  profiles  are  provided. 

Section  III  proposes  a  method  to  estimate  the  instantaneous  mass  of  the  nose  cone  that  is  suitable  for  use 
in  a  dynamic  simulation.  In  Section  IV,  a  method  is  described  that  allows  one  to  estimate  the  location  of  the 
OML  and  aerodynamic  coefficients  between  test  or  prediction  points,  based  on  the  instantaneous  mass.  In 
Section  V  we  describe  the  method  used  to  generate  the  aerodynamic  models  for  the  example  body  described 
in  Section  VI. 


II.  Generation  of  Outer  Mold  Lines  from  Test  Data 

In  order  to  generate  an  aerodynamic  model  of  an  ablated  vehicle,  an  approximation  for  the  outer  mold 
line  must  be  constructed.  The  shape  of  nose-cones  from  experimental  recession  measurements  are  often 
provided  in  the  form  of  upper  and  lower  curves  that  describe  the  surface  of  the  body  in  the  nominal  plane  of 
symmetry,  i.e  the  body  x-z  plane.  One  objective  of  this  section  is  to  develop  a  method  of  estimating  out-of¬ 
plane  grid  points  so  that  an  aerodynamics  prediction  or  CFD  code  can  be  used  to  estimate  the  aerodynamic 
characteristics  .  For  cases  where  only  2-D  profiles  are  obtained  from  recession  measurements,  OML  grid 
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points  outside  of  the  x-z  plane  must  be  estimated.  The  method  presented  here  was  developed  for  nose  cones 
that  are  axissymmetric  before  ablation  occurs. 

Here  we  propose  to  estimate  out-of-plane  grid  points  by  interpolating  between  two  surfaces  of  revolution 
about  the  x-body  axis.  Figure  1  shows  how  the  upper  and  lower  curves  defining  a  body  in  the  x-z  plane  can 
be  used  to  generate  an  interpolated  body  that  can  be  used  to  obtain  estimates  of  OML  coordinates  in  the 
y-z  plane. 

Figure  1  shows  a  side  and  sectional  view  of  a  nose  cone.  The  requirements  are  that  the  interpolated  OML 
must  be  coincident  with  the  upper  curve  when  0  =  n/ 2  and  with  the  lower  curve  when  0  =  —  7r/2,  where  0 
is  the  angle  (in  the  y-z  plane)  measured  from  the  y-axis  to  a  point  on  one  of  the  surfaces.  Furthermore,  the 
interpolated  crossections  are  taken  to  be  smooth  at  6  =  ±tt/2  and  symmetric  about  the  x-z  plane.  These 
boundary  conditions  are  motivated  by  the  fact  that  the  OML  of  the  unablated  vehicle  could  be  described 
by  a  surface  of  revolution.  For  lifting  reentry  vehicles,  test  data  shows  that  leeward  surfaces  ablate  less 
than  the  windward  surfaces.  Thus,  the  vehicle  OML  does  not  remain  symmetric  when  flying  at  a  non-zero 
angle-of- attack  under  thermal  conditions  where  ablation  occurs.  We  now  construct  an  interpolation  function 
that  satisfies  the  position  and  slope  boundary  conditions  described  above. 

Let  A (6)  E  [0, 1]  be  an  interpolating  function  such  that: 

ri(6)  =  \(0)r u  +  [1  -  X(6)]ri  (1) 

where  is  the  radius  of  the  interpolated  surface,  ru  is  the  radius  of  upper  surface  (leeward  side)  and  77  is 
radius  of  the  lower  surface  (windward  side) .  It  is  important  to  note  that  the  interpolating  function  is  only  to 
be  used  on  the  interval  6  E  [— 7r/2,7r/2]  and  that  construction  of  the  interpolated  surface  in  the  2nd  and  3rd 
quadrants  of  the  y-z  plane  must  be  made  using  symmetry  in  quadrants  1  and  4.  Note  that  a  simple  linear 
interpolating  function  will  not  satisfy  the  continuity  conditions  because  a  linear  function  does  not  contain 
enough  free  parameters  to  satisfy  both  the  position  and  slope  constraints  at  0  =  ±tt/2  .  A  cubic  polynomial 
interpolating  function;  however,  will  satisfy  the  four  boundary  conditions,  i.e. 

ri(  tt/2)  =  ru 

n(-7r/2)  =  n  (o\ 

^/2)  =  ° 

The  general  form  for  a  third  order  interpolating  polynomial  is  given  by: 

A(6>)  =  a03  +  b02  +  c0  +  d  (3) 


The  coefficients  of  the  polynomial  equation  must  satisfy  the  following  equation: 


'  (tt/2)3  (tt/2)2  (tt/2)  1' 

a 

"l" 

(-tt/2)3  (-tt/2)2  (-tt/2)  1 

b 

0 

3(tt/2)2  2(tt/2)1  1  0 

c 

0 

_3(— 7r/2)2  2(- tt/2)1  1  0_ 

_d_ 

_0_ 

Solving  for  the  polynomial  coefficients  one  obtains: 


a 

‘-2/tt3‘ 

b 

0 

c 

3/(27 r) 

A 

.  1/2  . 

(4) 


(5) 


Substituting  the  results  of  Equation  5  back  into  Equation  1  one  can  estimate  the  location  of  points  on 
the  outer  mold  line  in  the  y-z  plane  at  a  given  x-station  for  known  values  of  ru  and  77  via 


Vi  =  n  cos(6>) 

Zi  =  ri  sin(0)  (6) 
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Figure  1.  Cubic  interpolation  between  two  bodies  of  revolution 


3  of  18 


American  Institute  of  Aeronautics  and  Astronautics 


The  geometric  properties  of  such  an  interpolated  profile  will  be  of  interest  in  subsequent  sections.  In 
particular,  the  area  bounded  by  the  region  defined  by  Equation  1  with  polynomial  coefficients  defined  by 
Equation  5  is  given  by: 


A  =  J  ^  rj(0)dd  =  |r(13 rf  +  13 r2u  +  9 rurt)  (7) 

It  will  also  be  useful  to  compute  the  volume  of  regions  bounded  above  and  below  by  surfaces  parallel  to 
the  y-z  plane  enclosed  by  Equation  6  and  linearly  interpolated  in  x  as  shown  in  Figure  2.  Such  a  region  will 
be  used  as  an  estimate  of  the  shape  of  an  ablated  vehicle  between  any  two  x  body  axis  stations  for  which 
data  exists  in  the  form  of  rUl ,  rq ,  aq,  rU2,  n2,  x<i.  Without  loss  of  generality,  we  will  assume  that  the  base  of 
a  section  lies  at  x\  =  0.  We  can  then  obtain  an  expression  for  the  interpolated  radius  at  any  location  (#,  x) 
as 

n(0,x)  =  -T( n(e,X2 )  -  ri(9,0))(x)  +  ri(9, 0)  (8) 


Figure  2.  Volume  enclosed  by  a  section  of  interpolated  nose  cone 


The  volume  of  the  region  can  be  computed  as: 
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Substituting  the  interpolation  coefficients  of  Equation  5  into  Equations  1  and  8  and  evaluating  the  integral 
one  can  obtain: 


V=^- 

3 


Ax+  A2  +  —  [26(  ^U2^Ui  +  rhrh)  +  9(r  u2rh  +  n2r  uj] 


(10) 


where  A\  and  A 2  are  the  areas  of  the  base  and  top  of  the  solid  respectively. 

Finally,  the  reader  should  be  advised  that  there  are  limitations  of  the  cubic  radial  interpolation  function 
used  in  this  application.  In  order  for  the  interpolation  function  to  produce  a  physically  meaningful  result  at 
a  given  x-station,  the  lower  radius  must  be  no  less  than  about  38%  of  the  upper  radius  or  more  precisely, 


ri  > 


6 


7T2  +  6 


(ii) 


If  this  relationship  is  not  satisfied,  then  the  cubic  interpolation  function  will  result  in  cusped  cross-section. 
The  above  relationship  can  be  found  by  using  a  2nd  derivative  test  for  a  critical  point  at  0  =  —  7r/2,  i.e. 


d2z  _  d2z  cPO 
dy 2  ^tt  d02  dy 2 


=  0 


Defining  n  =  kru  and  substituting  Equation  1  into  Equation  12,  one  obtains: 


0  = 


d2z 

dO2 


d2  [X (0)ru  +  (1  -  X(0))]kru]  sin(0) 


dO2 


evaluating  Equation  13  and  solving  for  k  yields, 


k  = 


6 


7T2  +  6 


(12) 


(13) 


(14) 


□ 


For  cases  where  the  windward  vehicle  surface  ablates  significantly  more  than  the  leeward  surface,  the 
reader  is  advised  to  check  this  condition  at  each  x-station  to  ensure  that  the  interpolation  function  provides 
a  physically  meaningful  OML. 

III.  Computing  Instantaneous  Mass  of  Nose  Cone 

First  we  will  assume  that  the  nose  cone  density  pn  is  constant  throughout.  It  is  also  assumed  that  test 
data  is  available  that  provides  the  upper  and  lower  curves  for  one  or  more  ablated  cases.  The  volume  and 
thus  the  mass  of  the  ablated  and  unablated  nose  cones  can  computed  using  the  3-D  OMLs  generated  using 
the  techniques  described  in  the  previous  section.  The  following  formula  for  predicting  mass  loss  rates  due  to 
ablation  has  been  proposed:3 

dm  Ch  1  ta3  ( 


-  = —  ~pVc 

dt  H  2^ 


S 


(15) 


where  m  is  the  mass  of  the  nose  cone,  Ch  is  a  dimensionless  heat  transfer  coefficient  less  than  unity,  H  is 
the  effective  heat  of  ablation,  p  is  the  air  density,  V  is  the  velocity  of  the  body  with  respect  to  the  air  mass, 
and  S  is  the  current  cross  sectional  frontal  area  of  the  body.  Note  that  Ch  and  H  are  vehicle  and  material 
dependent  parameters  and  we  hereby  denote  the  ratio  Ch/H  as  Kh •  Equation  15  can  be  integrated  over 
trajectory  to  obtain: 


to/  -  m|  =  -\-KhS  f  p(t)V(t)3dt 

z  Jo 


(16) 


where  mj  and  are  the  initial  and  final  values  of  the  nosecone  mass.  Denoting  Am  =  mf  —  mi  and  solving 
for  Kh  we  have: 

Kh  =  2 - -  (17) 

S  Jo  p{t)V(t)3dt 

Estimates  of  p  can  be  obtained  from  an  atmospheric  model,  V  can  be  measured  or  estimated  from  test 
instrumentation,  and  changes  in  mass  can  be  determined  from  recession  measurements.  Thus,  Equation  17 
can  be  used  to  empircally  estimate  Kh  for  the  flight  profile  for  which  the  recession  measurements  were 
obtained.  The  integral  is  not  analytically  tractable;  however,  it  can  be  evaluated  by  quadrature. 
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IV.  Interpolating  Outer  Mold  Lines 


While  the  centroid  and  elements  of  the  inertia  tensor  of  the  nose  cone  can  change  significantly  as  a 
result  of  ablation,  those  effects  usually  have  little  influence  upon  the  mass  properties  of  the  vehicle  as  a 
whole.  In  some  cases  the  changes  in  mass  properties  due  to  ablation  can  be  neglected  for  the  purposes  of 
flight  simulation.  Nevertheless,  computation  of  the  OML  based  on  instantaneous  mass  is  still  important 
because  the  aerodynamic  properties  are  a  strong  function  of  the  OML.  Since  the  aerodynamic  properties 
are  primarily  influenced  by  OML,  and  OML  and  mass  are  related,  a  method  is  proposed  that  will  allow 
the  aerodynamic  characteristics  of  an  ablated  profile  to  be  estimated  using  using  instantaneous  mass  as  an 
interpolation  variable. 

Generally,  a  small  number  of  OMLs  are  known  from  test  or  prediction  data.  We  would  like  to  estimate 
the  shape  and  aerodynamic  characteristics  of  the  vehicle  based  on  measurements  or  predictions  of  a  finite 
number  of  OMLs  and  estimates  of  instantaneous  mass.  During  a  simulation  run,  the  time  rate  of  change  of 
vehicle  mass  due  to  ablation  can  be  computed  dm/dt  via  Equation  15.  Thus  the  mass  of  the  vehicle  can  be 
computed  as  a  function  of  time  in  a  dynamic  simulation.  Given  the  masses  associated  with  known  OMLs, 
it  is  desired  to  estimate  the  instantaneous  OML  location  and  the  associated  aerodynamic  characteristics  for 
any  specified  mass.  In  order  to  obtain  a  point  on  an  OML  for  a  given  mass,  one  must  interpolate  between 
two  points  on  adjacent  OMLs  defined  by  test  or  prediction  data.  Numerous  choices  exist  for  selecting  match 
points  on  the  inner  and  outer  surfaces.  Here  we  propose  to  match  points  on  adjacent  OMLs  that  have 
equal  slopes.  This  choice  is  motivated  by  Newtonian  flow  considerations  wherein  the  local  temperature  and 
pressure  on  a  body  are  dictated  by  the  local  surface  inclination  measured  with  respect  to  the  freestream,4 
i.e. 


T2 


=  27(7  - 1) 
'  (7+1)2 

A  =  JA_ 

Poo  7  +  1 


sin2  (3 
Ml  sin2  p 


(18) 


where  /?  is  the  angle  between  a  tangent  line  on  a  surface  and  the  freestream.  We  therefore  propose  to 
interpolate  along  lines  connecting  points  on  the  inner  and  outer  surfaces  that  have  equal  surface  inclination 
angles  as  illustrated  in  Figure  3.  Points  on  the  interpolated  OML  are  computed  as  follows: 


XOML  =  A m(m)xi  +  (1  -  A m{m))x0 
zqml  =  A m{m)zi  T  (1  Am (m) )zG 


(19) 


where  A m(ra)  G  [0, 1]  is  an  interpolation  function  that  relates  changes  in  distance  to  changes  in  the  instan¬ 
taneous  mass  of  the  vehicle.  In  cases  where  discontinuities  arise  because  segments  of  one  OML  do  not  have 
slopes  that  match  those  on  adjacent  surface,  such  points  are  matched  to  a  corner  point  on  an  adjacent  OML 
as  illustrated  in  Figure  3. 

Since  mass  or  volume  is  a  cubic  function  of  the  linear  dimensions  of  a  body,  one  might  postulate  the 
following  form  for  the  interpolation  function: 

A  m(ra)  =  K\mfc  +  K2  (20) 


Clearly  A m(raG)  =  0  and  A m{mi)  =  1?  thus 


Kx 


K2 


1 

i  i 

mf  —  ml 

-ml 
t  r 

mf  —  ml 


(21) 


The  distance  that  a  point  moves  along  a  line  connecting  the  inner  and  outer  control  points  is  therefore 
assumed  to  be  related  to  the  cube  root  of  the  volume  or  mass.  The  above  relationship  would  yield  accurate 
results  for  an  object  that  ablated  symmetrically  in  all  directions;  however,  ablation  on  reentry  vehicles 
is  typically  asymmetric.  One  should  therefore  judiciously  select  an  interpolation  function  that  correlates 
changes  in  mass  with  changes  in  OML  location.  For  example,  if  an  unablated  body  is  a  cylinder  capped  by 
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Figure  3.  OML  Interpolation  based  on  instantaneous  mass 


a  hemisphere  and  held  at  zero  angle  of  attack  in  a  hypersonic  flow,  one  would  expect  the  principle  change  in 
geometry  to  be  in  the  length  of  the  body  with  perhaps  a  slight  change  in  radius.  Since  the  mass  or  volume 
of  such  a  body  is  a  linear  function  of  length,  one  would  expect  that  a  linear  interpolation  function  would  be 
an  appropriate  choice  for  relating  changes  in  mass  to  OML  location. 


In  general, 


A  m(m)  =  Kim  +  K2 

1 

K,= 


rrii  -  mQ 
-m„ 


A  m(ra)  =  Kif{m)  +  K2 


(22) 

(23) 


(24) 


where  /(m)  is  a  function  that  correlates  mass  and  OML  location,  and  the  coefficients,  K\  and  K2  are 
computed  from  the  boundary  conditions  A m(ra0)  =  0  and  Am(rn^)  =  1. 

The  interpolation  of  aerodynamic  tables  using  instantaneous  mass  proceeds  along  similar  lines  since 
the  aerodynamic  properties  are  significantly  influenced  by  OML.  Typically  aerodynamic  force  and  moment 
coefficients  are  tabulated  as  a  function  of  angle-of-attack,  sideslip  angle,  Mach  number  and  control  surface 
deflections.  For  a  vehicle  that  experiences  ablation,  the  aerodynamic  data  will  also  be  a  function  of  the 
instantaneous  outer  mold  line  of  the  vehicle.  Thus,  additional  aerodynamic  tables  can  be  constructed  for 
a  representative  set  of  ablated  OMLs.  Between  the  table  break  points  established  for  given  or  estimated 
OMLs  it  is  proposed  that  the  tables  be  interpolated  using  the  interpolation  function  that  is  used  to  compute 
OML  location  based  on  mass,  i.e.  Equation  24,  because  of  the  strong  influence  of  OML  on  aerodynamic 
coefficients.  For  an  arbitrary  aerodynamic  coefficient  the  following  interpolation  scheme  is  proposed: 


CX(M,  <a,  /?,  <5,  m)  =  a,  (3,  S,  rrq)  +  [1  -  A (m)]Cx(M,a,/3,5,m0 


(25) 


In  summary,  to  interpolate  the  aerodynamic  tables,  one  must  estimate  the  instantaneous  mass  of  vehicle 
using  Equation  15,  and  use  an  interpolation  function  given  by  Equation  24  in  conjunction  with  Equation  25. 
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V.  Prediction  of  Aerodynamic  Properties  from  Outer  Mold  Lines 


Many  design  codes  exist  that  can  predict  the  aerodynamic  characteristics  of  blunted  bodies  at  high 
speeds.  Newtonian  theory  has  been  shown  to  provide  good  lift  and  pitching  moment  predictions  for  a  variety 
of  sharp  and  blunted  body  shapes  at  Mach  numbers  as  low  as  2. 75. 5  Addition  of  viscous  effects  (skin  friction, 
base  drag)  improves  the  accuracy  of  drag  predictions.  Newtonian  impact  theory  (Cp  =  2  sin2  S  )  typically 
provides  reasonable  agreement  for  sharp  cones.  The  modified  form  of  Newtonian  theory,  where  the  pressure 
coefficient  of  2  is  replaced  by  the  stagnation  pressure  behind  a  normal  shock,  has  been  shown  to  give  better 
results  for  highly  blunted  configurations. 

Although  Missile  Datcom6  contains  a  Newtonian  flow  option,  it  was  not  selected  for  the  analysis  of  the 
example  body  presented  in  the  subsequent  section,  because  it  is  limited  to  bodies  of  circular  or  elliptic  cross- 
section,  and  the  cross-section  shapes  resulting  from  the  proposed  ablateness  model  are  neither.  In  addition, 
the  Newtonian  formulation7  used  by  Missile  Datcom  does  not  include  the  axial  force  contribution  to  pitching 
moment  (a  common  simplification) .  The  code  selected  for  the  present  analysis  is  a  modified  version  of  the  Low 
Observable  Design  Synthesis  Tool8  (LODST),  which  was  developed  for  analysis  of  arbitrary  body  missiles. 
LODST  has  a  geometry  modeler  suitable  for  shapes  generated  by  the  proposed  ablateness  model,  and  also 
contains  a  Newtonian  flow  calculation.  Both  of  these  were  taken  from  the  Supersonic/Hypersonic  Arbitrary 
Body  Program,9  a  widely  used  code.  The  Newtonian  flow  calculations  were  isolated  from  the  slender  body 
theory  calculations  within  LODST,  and  base  drag  and  skin  friction  calculations  were  added  using  methods 
from  Missile  Datcom. 

Many  databases  exist  on  blunted  cones  at  high  speeds.  Neal10  tested  several  configurations  at  a  Mach 
number  of  6.77,  one  of  which  is  similar  to  the  configuration  considered  in  the  subsequent  example.  The  length 
to  base  diameter  ratio  of  the  configuration  considered  in  the  example  in  the  next  section  decreases  from  1.52 
to  0.98  as  the  nose  ablates.  In  order  to  validate  the  accuracy  of  the  predictions  for  the  example  problem, 
we  compare  results  obtained  from  the  LODST  code  to  test  data.  The  comparable  ratio  of  the  configuration 
tested  by  Neal  is  1.38.  The  ratio  of  the  nose  bluntness  radius  to  base  radius  of  the  example  configuration 
increases  from  0.44  to  0.73  as  the  nose  ablates.  The  comparable  ratio  of  the  configuration  tested  by  Neal  is 
0.66.  The  LODST  representation  of  the  Neal  configuration  is  shown  in  Figure  4.  A  comparison  of  the  axial 
force  coefficient  predicted  by  LODST  with  test  data  is  shown  in  Figure  5.  The  agreement  is  excellent  when 
the  viscous  contribution  is  added  to  the  modified  Newtonian  result.  Normal  force  comparisons  are  shown  in 
Figure  6.  Again,  the  agreement  is  excellent  for  the  modified  Newtonian  prediction. 

VI.  Example 

Consider  the  unablated  nose  cone  of  a  reentry  vehicle  shown  in  Figure  7  that  consists  of  a  12  deg  conical 
frustrum  capped  by  a  zone  of  a  sphere  of  radius  1.  More  precisely  the  outer  mold  lines  in  the  x-z  plane  are 
given  by: 


z(x)  =  ±\/l  —  x2i  cos(78°)  <  x  <  1  (26) 

z(x)  =  ±[tan(12°)x  —  (sin(78°)  +  tan(12°)  sin(12°))],  —  6  <  x  <  cos(78°)  (27) 

Suppose  that  at  a  certain  stage  in  a  flight,  a  recession  sensor  suite  has  provided  measurements  of  the  ablated 
nose  cone  in  the  x-z  plane  that  has  the  geometry  represented  by  inner  curves  shown  in  Figure  7  which  are 
described  as  follows: 


Zu(x) 

\J — 2(x  +  2), 

-  3.096  <  x  <  -2 

(28) 

-  = 

[—  tan(12°)x  +  (sin(78°)  +  tan(12°)  sin(12°))]  —  .2, 

—  6  <  x  <  —3.096 

(29) 

Zl{x )  = 

-3.714  <  x  <  -2 

(30) 

Zl(x)  = 

[tan(12°)x  —  (sin(78°)  +  tan(12°)  sin(12°))]  +  .3, 

-  6  <  x  <  -3.714 

(31) 

Based  on  the  outer-mold-lines  measured  in  the  x-z  plane  and  an  estimate  of  the  instantaneous  mass  of  the 
nose  cone  it  is  desired  to  estimate  the  instantaneous  aerodynamic  properties  and  OML  location.  This  is 
accomplished  by  using  Equations  19  and  25.  Points  on  the  ablated  and  unablated  surfaces  whose  tangents 
are  equal  are  used  in  Equation  19,  thus,  mass  interpolated  points  on  the  OML  in  the  x-z  plane  proceed  along 
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Figure  4.  LODST  Grid  for  Neal  Configuration 


the  equi-tangent  lines  shown  in  Figure  8.  Figure  9  shows  the  estimated  outer  mold  line  locations  associated 
with  specified  values  of  the  mass  interpolation  function  given  by  Equation  24. 

For  a  given  value  of  nose  cone  mass,  mass  interpolated  curves  in  the  x-z  plane  for  the  upper  and  lower 
surfaces  can  be  generated  as  described  above  and  one  can  use  Equations  1  and  6  to  estimate  out-of-plane  grid 
points  that  are  suitable  for  use  with  CFD  or  an  aerodynamic  prediction  code  such  as  LODST.  The  3-D  OML 
of  the  ablated  body  is  shown  in  Figure  10.  Upon  generating  such  a  grid,  the  aerodynamic  properties  can  be 
estimated  using  CFD  or  prediction  codes  such  as  LODST.  The  actual  volume  or  mass  of  the  interpolated 
body  can  be  calculated  by  discretizing  the  along  the  x-body  axis  and  summing  the  volumes  of  each  section 
given  by  Equation  10.  In  this  example,  we  assume  that  the  nose  cone  material  is  homogenous  with  unit 
density,  thus  the  magnitudes  of  the  volume  and  the  mass  are  equal.  Mass  calculations  can  be  used  to 
assess  the  adequacy  of  basis  functions  chosen  for  the  mass  interpolation.  The  effectiveness  of  the  outer-mold 
line  interpolation  method  can  be  assessed  by  specifying  a  value  for  the  interpolation  function  A,  estimating 
the  outer  mold  line  location  and  then  calculating  the  mass  of  the  body  described  by  the  resulting  outer 
mold  line.  The  mass  of  the  example  body,  which  we  will  call  the  actual  mass  (raa),  was  calculated  for  a 
number  of  values  of  A.  The  result  is  shown  in  Figure  11  and  represents  the  shape  of  the  ideal  interpolating 
function  for  the  geometry  under  consideration.  When  choosing  a  basis  function  for  Equation  24  one  can 
evaluate  the  effectiveness  of  the  choice  by  comparing  the  relationship  between  mass  and  A  for  a  given  form 
of  interpolation  function  to  the  ideal  result.  The  mass  associated  with  a  specified  value  of  A  can  be  obtained 
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Figure  5.  Axial  Force  Comparision  for  Neal  Configuration 


by  solving  Equation  24  for  m  which  we  will  define  as  ms : 


ms  =  f  1 


X-K2 

K1 


(32) 


From  Figure  11,  one  can  see  that  the  linear  and  cube-root  interpolation  methods  only  agree  with  the  ideal 
or  actual  values  at  the  boundary  conditions  (i.e.  unablated  and  fully  ablated  conditions)  and  that  the  linear 
interpolation  method  is  superior  to  the  cube-root  result  in  this  particular  case.  In  this  example,  we  note  that 
the  primary  change  in  mass  is  due  to  a  change  in  length  along  the  x-axis  and  that  there  is  relatively  little 
change  in  the  radial  direction.  Thus,  representing  the  change  in  outer  mold  line  location  as  a  linear  function 
of  mass,  yields  superior  results  to  using  cube-root  interpolation  since  the  major  changes  to  linear  dimensions 
of  the  body  occur  along  a  single  direction.  If  a  body  ablated  uniformly  in  all  directions,  one  would  expect 
the  cube-root  based  interpolator  to  produce  superior  results.  In  order  to  alleviate  computational  burden,  it 
is  desirable  to  avoid  running  CFD  or  aerodynamic  prediction  codes  at  every  time-step  in  a  3  or  6  degree  of 
freedom  dynamic  flight  simulation.  Thus  using  Equation  25  to  estimate  aerodynamic  properties  is  preferable 
to  generating  OML  grids  and  running  prediction  codes  as  the  simulation  progresses  even  though  it  is  expected 
that  the  latter  would  produce  more  accurate  results.  A  good  choice  of  interpolation  basis  function  can  reduce 
these  inaccuracies.  To  illustrate  this  point  we  use  LODST  to  predict  the  lift,  drag  and  pitching  moments  of 
the  ablated  and  unablated  bodies  at  Mach  =  10  and  a  =  10°.  Here  we  are  interested  in  assessing  how  well 
a  mass  interpolated  (i.e.  Equation  25)  aerodynamic  prediction  matches  the  LODST  prediction  for  OMLs 
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Figure  6.  Normal  Force  Comparision  for  Neal  Configuration 


associated  with  a  specified  mass.  For  the  purposes  of  this  exercise,  we  take  the  LODST  prediction  to  be  the 
”  actual”  value  of  the  aerodynamic  coefficient  which  we  use  to  assess  the  accuracy  of  Equation  25  for  both 
cube-root  and  linear  interpolation  methods.  Figures  12  and  ??  show  how  the  lift  and  drag  of  the  example 
body  change  as  a  function  of  mass.  Note  that  there  are  two  sources  of  error  associated  with  the  interpolated 
results  on  these  figures.  The  first  source  of  error  is  that  there  is  a  difference  between  the  interpolated  value 
of  mass  and  the  actual  mass  of  the  body  described  by  the  OML  grid,  i.e.  ms  —  ma  ^  0.  The  second  source 
of  error  arises  because  the  interpolation  basis  functions  do  not  yield  exact  results  at  points  other  than  the 
boundary  conditions.  This  example  shows  a  rather  extreme  case  where  there  is  a  large  difference  between 
two  known  OMLs.  In  practice  the  differences  between  known  OMLs  is  generally  much  smaller,  making 
the  choice  of  interpolation  basis  functions  less  critical.  Even  in  cases  where  large  differences  exist  between 
known  OMLs,  the  use  of  OML  interpolation  in  conjunction  with  prediction  codes  can  be  used  to  obtain  more 
accurate  estimates  of  aerodynamic  coefficients  and  mass  between  test  data  points  in  order  to  augment  table 
lookup  data  for  dynamic  simulation  purposes. 

VII.  Summary 

A  method  for  estimating  the  effects  of  nose  cone  ablation  on  the  aerodynamic  forces  and  moments  of  a 
reentry  vehicle  has  been  presented.  The  method  provides  a  means  of  generating  a  rough,  but  representative 
model  that  makes  use  of  interpolated  test  or  prediction  data  using  instantaneous  nose  cone  mass  as  an  inter- 
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Figure  7.  Example  Geometry  of  an  Ablated  and  Unablated  Nose  Cone 
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Figure  8.  Lines  connecting  points  of  equal  tangents  on  ablated  and  unablated  surfaces 


polation  variable.  Since  aerodynamic  properties  are  strongly  influenced  by  outer-mold-line,  an  interpolation 
function  was  proposed  that  can  be  used  to  approximate  the  location  of  points  on  the  outer-mold-line  as  well 
as  the  aerodynamic  characteristics  test  points.  While  the  changes  in  overall  vehicle  mass  properties  due  to 
ablation  are  typically  small,  one  can  make  use  of  the  interpolated  outer-mold-line  information  to  estimate 
mass  properties,  such  as  moments  of  inertia  or  center  of  gravity.  The  interpolated  aerodynamic  and  mass 
properties  provide  a  means  for  generating  data  for  dynamic  simulations  for  use  in  vehicle  proof-of-concept 
studies. 
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Figure  10.  Grid  of  points  on  ablated  OML  using  radial  interpolation. 
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Figure  11.  Mass  versus  interpolation  function  value 
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Figure  12.  Predicted  lift  coefficient  variations  due  to  mass 
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Figure  13.  Predicted  drag  coefficient  variations  due  to  mass 
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